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Introduction: 


Our  goal  is  to  develop  a  diagnostic  device  for  breast  cancer  based  on  fluorescence  and  diffuse 
reflectance  spectroscopy.  Our  first  objective  include  was  to  characterize  the  fluorescence  properties  of  normal, 
benign,  and  malignant  breast  tissue  in  patients  undergoing  breast  cancer  and  breast  reduction  surgeries.  Our 
second  objective  was  to  characterize  the  intrinsic  fluorescence  properties  of  tissue  constituents  in  order  to 
identify  the  biological  basis  for  spectroscopic  differences  in  tissue.  An  initial  ex  vivo  study  has  been  completed, 
and  a  fast,  portable  device  will  be  employed  to  measure  the  fluorescence  and  diffuse  reflectance  properties  of 
tissue,  in  vivo,  in  the  coming  year.  Results  are  correlated  with  the  gold  standard,  histopathology,  to  develop 
statistical  methods  to  determine  the  diagnostic  potential  of  this  technology.  Also,  we  would  like  to  note  that 
Changfang  Zhu  has  taken  over  as  PI  on  this  grant. 

Body: 

The  first  major  focus  of  work  over  the  past  12  months  has  been  to  characterize  tissue  samples  ex  vivo  in 
accordance  with  parts  b-d  of  Task  1.  18  patients  participated  in  this  study,  from  which  we  collected  and 
characterized  47  tissue  samples  for  their  fluorescence  and  diffuse  reflectance  properties.  A  statistical  analysis 
was  then  carried  out  to  identify  differences  in  the  spectra  of  malignant  and  non-malignant  breast  tissues  and 
determine  the  diagnostic  efficacy  of  these  techniques.  It  was  found  that  optical  spectroscopy  in  the  ultraviolet- 
visible  spectral  range  shows  distinct  differences  between  malignant  and  non-malignant  breast  tissues,  and  using 
a  statistical  method,  this  technique  can  discriminate  malignant  from  non-malignant  breast  tissues  with  overall 
accuracy  of  83%,  sensitivity  of  60%  and  specificity  of  90%.  In  addition,  we  explored  the  effect  of  probe 
geometry  on  the  optical  spectroscopic  diagnosis  of  breast  cancer.  It  has  been  shown  that  the  attenuation  of 
spectral  intensities,  which  is  characterized  by  a  multi-separation  fiber  optic  probe  used  in  this  study,  can  be 
utilized  for  breast  cancer  diagnosis.  A  manuscript  detailing  this  work  has  been  prepared  and  is  attached  in 
appendix  A. 

The  second  focus  of  work  over  the  past  year  has  been  to  develop  a  system  capable  of  characterizing 

tissue  in  vivo.  To  this  end,  a  novel  side  firing  probe  has  been  designed  that  is  capable  of  making  measurement 
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through*  needles  currently  used  for  breast  core  needle  biopsy.  A  first  generation  side  firing  probe  has  been 
fabricated  and  experimentally  assessed  for  its  performance.  In  addition,  we  have  obtained  local  IRB  approval 
and  expect  to  collect  in  vivo  measurements  on  approximately  50  samples  in  the  coming  year,  which  will  satisfy 
part  b  of  Aim  2.  One  minor  change  we  would  like  to  make  is  that  we  have  decided  to  focus  on  making  in  vivo 
measurements  during  core  needle  biopsy  procedure,  rather  than  during  breast  cancer  surgeries.  Due  to  concerns 
about  extending  the  time  of  surgery  and  the  logistics  of  making  measurements  in  the  operating  room,  we  have 
concluded  that  it  is  much  more  feasible  for  us  to  make  in  vivo  measurements  on  patients  undergoing  core  needle 
biopsies,  rather  than  the  surgical  patients.  We  are  currently  in  the  process  of  trying  to  have  this  changed  in  our 
official  statement  of  work.  Much  of  the  work  for  Aim  3  has  also  been  accomplished  with  the  development  of 
analysis  techniques  employed  in  the  ex  vivo  study  and  described  in  detail  in  Appendix  A,  so  we  expect  to  fully 
complete  all  aims  by  the  end  of  our  funding  period. 

Key  Research  Accomplishments 

•  Identification  of  spectral  differences  present  in  human  breast  tissue,  ex  vivo. 

•  Development  of  a  statistical  method  for  analyzing  the  spectral  data  and  evaluating  the  classification 
accuracy  of  optical  spectroscopic  diagnosis  of  breast  cancer 

•  Development  of  a  system  capable  of  making  rapid  measurements  in  vivo,  during  core  needle  biopsy 
procedures 

Reportable  Outcomes: 

Journal  Articles: 

•  Changfang  Zhu,  Gregory  M.  Palmer.  Tara  M.  Breslin,  Fushen  Xu,  Nimmi  Ramanujam,  The  Use  of  a 
Multi-separation  Fiber-Optic  Probe  for  the  Optical  Diagnosis  of  Breast  Cancer,  submitted  to  Journal  of 
Biomedical  Optics  (2004) 

Conference  Presentations: 

•  Palmer,  G.M.,  C.  Zhu,  C.  Lubawy,  and  N.  Ramanujam.  “Exploiting  intrinsic  sources  of  optical  contrast 
for  cancer  detection.”  Oral  Presentation,  American  Chemical  Society  Annual  Meeting  and  Exposition, 
Philadelphia,  PA,  Aug.  2004. 
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,  •  Khu,  C.,  G.M.  Palmer,  T.M.  Breslin,  F.  Xu,  N.  Ramanujam,  “Fluorescence  Spectroscopic  Diagnosis  of 

Breast  Cancer  Using  a  Multi-separation  Fiber  Optic  Probe”,  Poster  Presentation,  2004  Gordon  Research 
Conference:  Lasers  in  Medicine  &  Biology,  Meriden,  NH,  Jul.  2004 

•  Palmer,  G.M.,  C.  Zhu,  M.C.  Skala,  T.M.  Breslin,  F.  Xu,  K.W.  Gilchrist,  Q.  Liu,  and  N.  Ramanujam.  “A 
model  of  diffuse  reflectance  and  its  application  to  breast  cancer  diagnosis.”  Oral  Presentation,  AAMI 
Annual  Meeting  -  Young  Investigator  Competition,  Boston,  MA,  June,  2004 

•  Palmer,  G.M.,  Q.  Liu,  N.  Ramanujam,  C.  Zhu,  T.M.  Breslin,  F.  Xu,  and  K.W.  Gilchrist.  “Monte  Carlo- 
based  inverse  model  of  diffuse  reflectance  for  determination  of  UV-VIS  optical  properties  and  its 
application  to  breast  cancer  diagnosis.”  Oral  Presentation,  Optical  Society  of  America,  Biomedical 
Optics  Topical  Meeting,  Miami,  FL,  April,  2004 

•  Palmer,  G.M.,  C.  Zhu,  T.M.  Breslin,  F.  Xu,  K.W.  Gilchrist,  Q.  Liu,  and  N.  Ramanujam.  “Physically  vs. 
empirically  based  algorithms  for  diagnosis  of  breast  cancer  using  optical  spectroscopy.”  Oral 
Presentation,  Biomedical  Engineering  Society  Annual  Meeting,  Nashville,  TN,  Oct.  2003. 

Conclusions 

Significant  differences  in  the  fluorescence  and  diffuse  reflectance  properties  of  normal  and  malignant  breast 
tissues  have  been  identified,  and  a  method  for  classification  of  these  tissues  has  been  developed.  Furthermore 
we  have  completed  design  of  a  fiber  optic  probe  capable  of  measuring  fluorescence  and  diffuse  reflectance  in 
vivo  during  breast  core  needle  biopsy,  and  expect  to  begin  this  in  vivo  study  soon.  A  diagnostic  device  based  on 
fluorescence  and  diffuse  reflectance  spectroscopy  has  the  advantage  of  being  fast,  quantitative,  and  minimally 
invasive,  and  has  the  potential  to  improve  breast  cancer  care  both  in  assessing  surgical  margins  for  residual 
cancer  in  vivo  during  surgery,  as  well  as  during  core  needle  biopsy  to  accurately  identify  the  appropriate  region 
for  biopsy. 
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Abstract 


This  study  explored  the  effects  of  the  illumination  and  collection  geometry  on  optical  spectroscopic 
diagnosis  of  breast  cancer.  Fluorescence  and  diffuse  reflectance  spectroscopy  in  the  ultraviolet-visible 
wavelength  range  were  made  with  a  multi-separation  probe  at  three  illumination-collection  separations  of  735, 
980  and  1225  pm,  respectively,  from  13  malignant  and  34  non-malignant  breast  tissues.  Statistical  analysis  was 
carried  out  on  two  types  of  data  inputs:  (a)  the  fluorescence  and  diffuse  reflectance  spectra  recorded  at  each  of 
the  three  illumination-collection  separations,  and  (b)  the  integrated  fluorescence  (at  each  excitation  wavelength) 
or  diffuse  reflectance  over  the  entire  spectrum  at  all  three  illumination-collection  separations.  The  results  from 
this  study  show  that  using  the  integrated  fluorescence  intensities  recorded  at  a  single  excitation  wavelength  at 
all  three  illumination-collection  separations  can  discriminate  malignant  from  non-malignant  breast  tissues  with 
similar  classification  accuracy  to  that  using  spectral  data  measured  at  several  excitation  wavelengths  with  a 
single  illumination-collection  separation.  These  findings  have  significant  implications  with  respect  to  the  design 
of  an  optical  system  for  breast  cancer  diagnosis.  Examining  the  intensity  attenuation  at  a  single  wavelength 
rather  than  spectral  intensities  at  multiple  wavelengths  can  significantly  reduce  the  measurement  and  data 
processing  time  in  a  clinical  setting  as  well  as  the  cost  and  complexity  of  the  optical  system. 

Keywords:  breast  cancer;  fluorescence;  reflectance;  spectroscopy;  probe  geometry 
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.1.  Introduction 


Surgical  biopsy  and  core  needle  biopsy  are  commonly  used  for  the  diagnosis  of  breast  lesions.  Compared  to 
surgical  biopsy,  core  needle  biopsy  is  less  invasive,  less  expensive,  faster,  and  requires  a  shorter  time  for 
recovery.  However,  its  sampling  accuracy  is  limited  because  only  a  few  small  pieces  of  tissue  are  sampled  from 
random  locations  in  the  suspicious  mass.  Consequently  the  needle  biopsy  procedure  has  a  false-negative  rate  of 
1%  -  7%  1  when  verified  with  follow  up  mammography,  and  repeat  biopsies  are  required  in  9%  -  18%  of 
patients  2’ 3. 

Optical  spectroscopy  has  been  increasingly  investigated  as  a  tool  for  breast  cancer  diagnosis  4.  In  particular, 
several  studies  have  demonstrated  that  there  are  significant  differences  in  the  ultraviolet-visible  (UV-VIS) 
fluorescence  and  diffuse  reflectance  spectra  of  normal,  benign  and  malignant  breast  tissues  5'14.  In  addition,  this 
technique  is  fast,  quantitative,  and  non-destructive.  There  are  important  benefits  to  be  reaped  by  incorporating 
optical  spectroscopy  as  an  adjunct  diagnostic  modality  to  core  needle  biopsy.  This  technology  can  be  deployed 
through  fiber-optic  probes  to  quickly  and  non-destructively  identify  the  tissue  type  (normal,  benign  and 
malignant)  at  the  needle  tip  during  a  breast  biopsy  procedure.  A  positive  reading  from  the  optical  measurement 
will  potentially  increase  the  likelihood  that  a  biopsy  is  being  sampled  from  a  tumor  site.  If  the  optical 
measurement  reads  negative,  then  the  needle  can  be  repositioned  (along  the  needle  track)  to  a  new  tissue  site. 
Currently  6-24  biopsies  are  taken  during  a  core  needle  biopsy  procedure.  If  the  optical  method  can  maximize 
sampling  from  tissue  sites  that  are  most  likely  to  be  cancerous,  and  minimize  unnecessary  removal  of  many 
normal  tissues,  it  could  make  the  breast  biopsy  procedure  more  accurate,  less  traumatic  to  the  patient  and  also 
reduce  the  number  of  biopsies  that  need  to  be  processed  in  order  to  obtain  a  confirmatory  diagnosis. 
Additionally,  if  optical  spectroscopy  proves  to  be  an  effective  way  of  identifying  cancerous  sites,  it  can  be 
incorporated  into  much  smaller  needles  than  the  11  gauge  Mammotome  needle  (for  example,  the  21  gauge 
needle  used  for  fine  needle  aspiration)  and  thus  make  this  procedure  as  minimally  invasive  as  possible. 
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.  Gupita  12  and  Majumder  et  al.  11  have  analyzed  different  spectral  data  sets  collected  from  the  same  set  of 
breast  tissues  ex  vivo  and  showed  that  the  fluorescence  emission  spectra  at  excitation  wavelengths  of  340  and 
488  nm  and  excitation  spectra  at  emission  wavelengths  of  390  and  460  nm  exhibit  differences  between  normal, 
benign  and  malignant  tissues.  Using  integrated  emission  intensities  at  340  nm  excitation  in  a  binary 
classification  scheme,  they  were  able  to  differentiate  malignant  from  normal  and  benign  tissues  with  a 
sensitivity  and  specificity  of  98%.  The  fluorescence  was  attributed  to  the  endogenous  fluorophores,  reduced 
nicotinamide  adenine  dinucleotide  (NADH)  and  collagen.  Yang  et  al.  5’ 7’ 8  showed  that  fluorescence  emission 
spectra  at  300  nm  excitation  and  excitation  spectra  at  340  nm  emission  could  be  used  to  discriminate  between 
malignant  and  normal  fibrous  tissues,  but  not  between  malignant  and  normal  adipose  tissues.  They  found,  for 
example,  that  the  ratio  of  normalized  intensities  at  268  and  289  nm  emission  (of  the  300  nm  excitation  spectra) 
discriminated  between  malignant  and  normal  fibrous  tissues  with  a  sensitivity  of  93%  and  specificity  of  95%. 
The  primary  endogenous  fluorophore  in  tissue  probed  at  these  excitation  and  emission  wavelengths  is 
tryptophan. 

Non-fluorescent  absorbers  and  scatterers  in  breast  tissues  also  contribute  in  part  to  differences  observed  in 
the  fluorescence  spectra  of  normal,  benign  and  malignant  tissues.  Diffuse  reflectance  provides  a  measure  of 
tissue  absorption  and  scattering.  Several  groups  have  explored  the  utility  of  diffuse  reflectance  spectroscopy 
between  250  -  800  nm,  for  breast  cancer  detection  ex  vivo  6’  8’  9’ 13  and  in  vivo  9’ 14.  Bigio  et  al.  measured  the 
diffuse  reflectance  spectra  through  a  core  biopsy  needle  and  during  breast  cancer  surgery  and  showed  that  this 
technique  can  differentiate  malignant  from  normal  tissues  with  a  sensitivity  of  60  -  70%  and  a  specificity  of  85 
-  95%  l4.  This  collection  of  studies  shows  that  diffuse  reflectance  spectra  can  be  used  to  differentiate  malignant 
from  normal  tissues.  Changes  in  diffuse  reflectance  of  malignant  tissues  are  likely  due  to  increased  protein  and 
hemoglobin  absorption,  increased  scattering  and  decreased  (3-carotene  absorption. 

In  these  previous  studies,  either  fluorescence  or  diffuse  reflectance  was  used  for  identification  of 
malignancy.  In  the  case  of  fluorescence,  spectra  were  measured  only  at  one  or  several  excitation  wavelengths. 
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Palmer  et  al.  10  measured  the  fluorescence  spectra  (at  a  total  of  nine  excitation  wavelengths  in  the  UV-VIS 
spectrum)  and  the  UV-VIS  diffuse  reflectance  spectra  of  freshly  excised  breast  tissues.  Using  a  multivariate 
statistical  algorithm,  they  found  that  four  out  of  the  ten  measured  spectra  were  sufficient  to  maximize  the 
discrimination  between  malignant  and  non-malignant  tissues.  These  included  emission  spectra  at  excitation 
wavelengths  of  300,  400,  420,  and  460  nm. 

The  study  described  in  this  paper  builds  upon  the  work  by  Palmer  et  al.  10.  The  primary  goal  of  this  study  is 
to  explore  the  effects  of  the  illumination  and  collection  geometry  on  the  optical  spectroscopic  diagnosis  of 
breast  cancer.  The  illumination  and  collection  geometry  is  an  important  component  of  tissue  optical 
spectroscopy  and  currently,  fiber-optic  probes  are  most  commonly  used  for  this  purpose.  In  previous  studies, 
fiber-optic  probes  with  only  one  fixed  separation  between  the  illumination  and  collection  fibers  were  employed 
for  optical  spectroscopy  of  the  breast 5"14.  These  probes  provided  a  fixed  optical  sensing  depth  in  tissue.  In  this 
study,  fluorescence  spectra  at  a  total  of  eight  excitation  wavelengths  and  UV-VIS  diffuse  reflectance  spectra 
were  measured  from  breast  tissues  with  a  multi-separation  fiber-optic  probe.  This  probe  has  three  illumination 
and  collection  separations,  which  enable  optical  spectra  to  be  measured  from  different  depths  within  the  tissue. 
The  fluorescence/reflectance  signal  attenuates  with  increasing  probing  depth  in  tissue  (due  to  the  increased  path 
length  of  the  photons).  Thus,  the  multi-separation  fiber-optic  probe  can  measure  the  attenuation  characteristics 
of  the  fluorescence  and  reflectance  measured  from  the  breast.  Two  different  types  of  data  sets  were  collected  in 
this  study  for  statistical  analysis:  (a)  the  fluorescence  and  diffuse  reflectance  spectra  recorded  at  each  of  the 
three  illumination-collection  separations  (or  synonymously,  different  optical  sensing  depths),  and  (b)  the 
fluorescence  (at  each  excitation  wavelength)  or  diffuse  reflectance  integrated  over  the  entire  spectrum  as  a 
function  of  illumination-collection  separation  (which  reflects  the  attenuation  of  light  within  the  tissue).  The 
statistical  analysis  of  both  types  of  data  sets  showed  that  the  integrated  fluorescence  intensities  at  a  single 
excitation  wavelength  recorded  at  all  three  illumination-collection  separations  discriminates  malignant  from 
non-malignant  breast  tissues  with  a  similar  classification  accuracy  compared  to  spectral  data  measured  at 
several  excitation  wavelengths  with  a  single  illumination-collection  separation.  These  findings  have  significant 
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implications.  Examining  the  intensity  attenuation  at  a  single  wavelength  rather  than  spectral  intensities  at 
multiple  wavelengths  (entire  spectrum)  could  significantly  reduce  the  acquisition  time  as  well  as  the  cost  and 
complexity  of  an  optical  system  for  breast  cancer  diagnosis.  This  study,  together  with  our  previous  study  10 
provides  important  information  on  the  utility  of  optical  spectroscopy  for  breast  cancer  diagnosis  during  core 
needle  biopsy. 

2.  Methods 

2.1.  Ex  vivo  sample  collection  and  handling 

The  breast  tissue  optical  spectroscopy  study  was  approved  by  the  Institutional  Review  Boards  at  the 
University  of  Wisconsin  -  Madison.  Breast  tissue  samples  were  obtained  from  patients  undergoing  either  a 
lumpectomy,  mastectomy  or  breast  reduction  surgery.  Fluorescence  and  diffuse  reflectance  spectra  were 
measured  on  the  freshly  excised  breast  tissues  within  2  hours  after  surgical  excision.  A  previous  study  carried 
out  by  our  group  15  showed  that  ex  vivo  spectra  measured  within  two  hours  after  excision  provides  a  relatively 
close  approximation  of  in  vivo  spectra.  During  the  tissue  measurement,  the  fiber-optic  probe  was  placed  in 
direct  contact  with  the  tissue  surface.  After  each  measurement,  the  exact  site  on  each  tissue  sample  where  the 
probe  was  placed  was  inked,  and  transverse  sections  were  cut  and  stained  for  histopathology.  Microscopic 
evaluation  was  performed  on  each  histological  section  by  a  board  certified  pathologist  (FX)  and  a  diagnosis  was 
established.  In  addition,  the  thickness  of  each  sample  was  measured  from  the  top  (where  the  ink  spot  was 
located)  to  the  bottom  of  the  H&E  stained  section.  Based  on  the  histological  breakdown,  each  sample  was 
broadly  classified  as  normal,  benign  or  malignant.  A  total  of  47  tissue  samples,  including  13  malignant,  32 
normal  and  2  benign  specimens,  were  collected  from  a  total  of  18  patients.  Table  1  shows  the  histological 
breakdown  of  the  47  samples  examined  in  the  breast  tissue  optical  spectroscopy  study.  For  samples  obtained 
from  lumpectomies  and  mastectomies,  the  average  tissue  thickness  was  5.6  ±  2.0  mm.  Sixteen  percent  of  the 
samples  had  a  thickness  of  ~  4  mm,  8%  had  a  thickness  of  ~  4.5  mm,  and  the  rest  had  a  thickness  greater  than  5 
mm.  All  samples  obtained  from  breast  reduction  surgeries  had  a  thickness  of  at  least  a  centimeter. 
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(Table  1) 


* 

2.2.  Optical  spectrometer  and fiber-optic  probe 

Fluorescence  and  diffuse  reflectance  spectra  were  measured  using  a  fiber-optic  probe  coupled  to  a  multi¬ 
wavelength  optical  spectrometer.  Figure  1  shows  the  (a)  schematic  of  the  optical  spectrometer  and  (b)  common 
end  of  the  fiber-optic  probe  that  comes  in  contact  with  the  tissue  (the  black  circles  correspond  to  the 
illumination  fibers,  the  gray  circles  correspond  to  the  collection  fibers  which  form  the  three  collection  rings,  and 
the  white  circles  correspond  to  dead  fibers  for  bundle  packing).  The  spectrometer  (Fig.  1(a))  consists  of  a  450 
W  xenon  lamp  (FL-1039,  J.Y.  Horiba),  a  scanning  double  excitation  monochromator  (Gemini  180,  J.Y. 
Horiba),  a  bifurcated  fiber-optic  probe  (Multimode  fibers),  a  filter  wheel,  an  imaging  spectrograph  (Triax  320, 
J.Y.  Horiba),  and  a  CCD  camera  (CCD3000,  J.Y.  Horiba).  The  common  end  of  the  fiber-optic  probe  (Fig.  1(b)), 
which  has  an  outer  diameter  of  3  mm,  consists  of  an  11 80  pm  diameter  illumination  core  and  three  concentric 
collection  rings  surrounding  the  core.  The  illumination  core  is  made  up  of  19  fibers,  each  of  which  has  a 
core/cladding  diameter  of  200/245  pm  and  a  numerical  aperture  (NA)  of  0.22.  Each  collection  ring  has  12  live 
fibers,  each  with  a  core/cladding  diameter  of  200/245  pm  and  an  NA  of  0.22.  The  remaining  fibers  are  dead 
fibers  of  the  same  size  for  bundle  packing.  The  illumination  diameter  maximizes  the  coupling  efficiency  from 
the  light  source,  and  the  signal-to-noise  ratio  for  the  fluorescence  measurements.  The  output  signals  from  the 
three  concentric  rings  of  collection  fibers  are  spatially  separated  on  the  CCD  chip,  thereby  allowing  for 
fluorescence  and  diffuse  reflectance  spectra  to  be  measured  at  three  illumination-collection  separations 
simultaneously.  The  illumination-collection  separation  is  defined  as  the  distance  from  the  center  of  illumination 
core  to  the  center  of  collection  fiber  in  each  ring.  The  three  illumination-collection  separations  in  this  fiber-optic 
probe  are  735,  980  and  1225  pm. 


(Figure  1) 
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.  Next,  Monte  Carlo  simulations  16  of  diffuse  reflectance  were  carried  out  to  verify  that  the  average  thickness 
of  tissue  samples  collected  from  patients  undergoing  lumpectomies  and  mastectomies  provides  a  semi-infinite 
geometry  for  diffuse  reflectance  spectroscopy  measurements  with  the  multi-separation  probe.  Simulations  were 
carried  out  on  a  homogeneous  tissue  model  with  absorption  coefficients  of  1  - 10  cm"1,  scattering  coefficients  of 
108  -  176  cm"1  and  an  anisotropy  factor  of  0.926.  The  range  of  optical  properties  used  is  within  that 
representative  of  breast  tissues  in  the  UV-VIS  spectrum  17.  The  simulated  medium  was  cylindrically  shaped 
with  a  thickness  of  4  mm,  which  represented  the  minimum  thickness  of  the  tissue  samples.  In  each  simulation, 
the  photons  were  launched  at  random,  uniformly  distributed  locations  over  a  range  of  angles  defined  by  an  NA 
of  0.22  and  over  a  circular  illumination  area  defined  by  the  fiber  bundle  diameter.  The  diffuse  reflectance 
escaping  the  medium  was  collected  over  a  circular  area  defined  by  the  fiber  diameter  and  over  a  range  of  exit 
angles  defined  by  an  NA  of  0.22.  The  refractive  index  above  the  medium  was  set  to  1 .452  to  simulate  an  optical 
fiber  and  that  below  the  medium  was  set  to  1.0.  The  refractive  index  of  the  medium  was  set  to  1.37.  The 
simulation  results  showed  that  the  light  transmittance  through  the  4  mm  thick  sample  was  less  than  5%  for  the 
range  of  optical  properties  used.  Specifically,  the  light  transmittance  through  the  sample  was  0.2%  when  the 
absorption  and  scattering  coefficient  were  10  and  176  cm"1,  respectively  and  4.8%  when  the  absorption  and 
scattering  coefficients  were  1  and  108  cm"1,  respectively.  These  results  indicate  that  the  average  tissue  thickness 
reported  in  this  study  provides  a  semi-infinite  geometry  for  diffuse  reflectance  spectroscopy  over  the  UV-VIS 
range.  The  fluorescence  emission  spectra  measured  in  this  study  occur  within  the  same  wavelength  range. 

To  assess  the  probing  depth  and  attenuation  of  the  fluorescence  measured  with  the  multi-separation  probe, 
Monte  Carlo  simulations  16  of  fluorescence  were  carried  out  on  a  homogeneous  tissue  model  with  a  wide  range 
of  optical  properties  in  the  UV-VIS  spectrum  (see  Table  2),  an  anisotropy  factor  of  0.926  at  both  the  excitation 
and  emission  wavelengths,  and  a  quantum  yield  of  unity.  The  simulated  medium  was  cylindrically  shaped  with 
a  thickness  of  5  mm.  It  should  be  pointed  out  that  in  real  tissue,  the  absorbers  and  fluorophores  have  distinct 
absorption  coefficients  and  the  fluorescence  efficiency  is  solely  related  to  the  product  of  the  absorption 
coefficient  and  quantum  yield  of  the  fluorophore.  However  in  the  Monte  Carlo  simulations,  these  two 
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.absorption  coefficients  are  combined  into  a  single  absorption  coefficient  and  the  fluorescence  efficiency  is  a 
function  of  the  overall  absorption  coefficient  (contributed  by  both  absorbers  and  fluorophores).  Thus  in  our 
Monte  Carlo  simulations,  when  the  absorption  coefficient  was  varied  at  the  excitation  wavelength,  the  quantum 
yield  was  appropriately  scaled  in  order  to  account  for  the  mismatch  in  fluorescence  efficiency  between  real 
tissue  and  the  simulated  medium  16. 

Table  2  shows  the  80%  probing  depth  and  the  percent  fluorescence  detected  with  the  three  illumination- 
collection  pairs  of  the  multi-separation  probe  from  a  homogeneous  fluorescent  medium,  with  (a)  variable 
absorption  coefficients  (pa)  and  a  fixed  scattering  coefficient  (ps),  and  (b)  variable  scattering  coefficients  and  a 
fixed  absorption  coefficient.  The  80%  probing  depth  is  defined  as  the  depth  above  which  80%  of  the  detected 
fluorescent  photons  originate.  The  percent  fluorescence  detected,  which  is  a  measure  of  the  collection 
efficiency,  is  the  number  of  fluorescent  photons  collected  with  a  single  fiber  in  each  collection  ring,  divided  by 
the  total  number  of  incident  photons  (5  million).  In  Table  2,  the  80%  probing  depth  varies  from  400  pm  to  1500 
pm.  For  a  medium  with  fixed  optical  properties,  the  probing  depth  increases,  while  the  percent  detected 
fluorescence  decreases  with  increasing  illumination-collection  separation.  For  a  fixed  illumination-collection 
separation,  the  probing  depth  decreases  with  increasing  absorption  or  scattering,  while  the  percent  detected 
fluorescence  decreases  with  increasing  absorption  and  increases  with  increasing  scattering,  respectively. 

In  summary,  for  a  given  set  of  optical  properties,  the  probing  depth  increases  only  modestly  with  increasing 
illumination-collection  separation,  while  the  percent  fluorescence  decreases  significantly  with  increasing 
illumination-collection  separation,  particularly  for  higher  absorption  and  scattering  coefficients.  Thus,  the  multi¬ 
separation  probe  geometry  can  be  used  to  measure  the  attenuation  characteristics  of  the  fluorescence  from 
breast  tissues.  It  is  expected  that  this  probe  geometry  will  measure  similar  trends  in  diffuse  reflectance  of  the 
breast.  One  of  the  objectives  in  this  paper  is  to  explore  whether  the  fluorescence/reflectance  attenuation  with 
increasing  probing  depth  can  be  utilized  for  discriminating  malignant  from  non-malignant  breast  tissues. 
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(Table  2) 


2.3.  Optical  spectroscopy  of  breast  tissues 

Fluorescence  emission  spectra  were  recorded  at  8  excitation  wavelengths  ranging  from  300  to  440  nm,  in  20 
nm  increments.  The  slit  width  of  the  excitation  monochromator  (Gemini  180)  was  set  to  provide  an  excitation 
band  pass  of  6.2  nm.  The  output  power  at  the  common  end  of  the  fiber-optic  probe  was  measured  using  a  low 
power  detector  (Newport  818-UV,  Newport)  connected  to  a  hand-held  optical  power  meter  (Newport  840-C, 
Newport)  and  ranged  from  10-20  pW  over  the  300  to  440  nm  excitation  wavelength  range.  At  each  excitation 
wavelength,  fluorescence  emission  spectra  were  measured  over  a  260  nm  wavelength  range,  with  the  first 
wavelength  shifted  by  20  nm  from  the  excitation  wavelength.  The  slit  width  of  the  imaging  spectrograph  was 
set  to  provide  a  spectral  resolution  of  7.9  nm.  Each  intensity-wavelength  point  in  the  emission  spectrum  was 
binned  over  1 8  pixels  on  the  CCD  chip,  resulting  in  a  wavelength  increment  of  4.7  nm.  The  integration  time  for 
each  fluorescence  spectrum  ranged  from  0.12  s  to  2  s,  depending  on  the  excitation  wavelength  used.  The  diffuse 
reflectance  spectra  were  recorded  from  350  to  600  nm.  The  slit  widths  of  the  excitation  monochromator  and 
imaging  spectrograph  were  chosen  to  provide  a  band  pass  of  3.5  nm  and  7.9  nm,  respectively.  Each  intensity- 
wavelength  point  in  the  diffuse  reflectance  was  binned  over  18  pixels  on  the  CCD  chip,  resulting  in  a 
wavelength  increment  of  4.7  nm.  The  integration  time  for  each  diffuse  reflectance  spectrum  was  0.05  s.  The 
overall  acquisition  time  was  approximately  1  minute  for  a  single  scan  of  8  fluorescence  spectra  and  one  diffuse 
reflectance  spectrum. 

All  fluorescence  and  diffuse  reflectance  spectra  were  calibrated  in  order  to  correct  for  the  (1)  background 
spectrum,  (2)  wavelength  dependence,  and  (3)  throughput  of  the  system.  The  background  spectrum,  which  was 
measured  with  the  probe  immersed  in  distilled  water  using  the  same  experimental  setup  for  optical  spectroscopy 
of  tissue,  was  first  subtracted  point-by-point  from  each  spectrum  prior  to  further  calibration.  The  fluorescence 
emission  spectral  intensities  at  each  excitation  wavelength  were  normalized  to  the  output  excitation  power,  to 
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qccount .for  the  wavelength-dependent  variation  of  the  excitation  light  intensity.  The  fluorescence  spectrum  was 
then  corrected  for  the  wavelength  dependent  response  of  the  collection  fibers,  imaging  spectrograph  and  CCD 
camera,  by  multiplying  it  point-by-point  by  the  correction  factors  measured  using  a  NIST  tungsten  calibration 
lamp.  This  procedure  also  corrected  for  differences  in  the  collection  efficiency  of  each  collection  ring.  Finally, 
each  fluorescence  spectrum  was  divided  by  the  peak  fluorescence  intensity  (excitation-emission  wavelength  of 
460-580  nm)  measured  with  the  probe  placed  on  the  face  of  a  quartz  cuvette  containing  a  solution  of 
Rhodamine  B  (2.14pM)  dissolved  in  ethylene  glycol  (115H3423,  Sigma  Chemical  Co.,  MO)  using  the  same 
integration  time,  to  account  for  the  time-dependent  changes  in  the  throughput  of  the  instrument.  The  diffuse 
reflectance  spectrum  was  calibrated  for  the  wavelength-dependent  response  and  the  throughput  of  the  system  by 
normalizing  it  to  the  diffuse  reflectance  spectrum  measured  with  the  common  end  of  the  fiber-optic  probe 
inserted  into  an  integration  sphere  (DRA-CA-30I,  Labsphere,  Inc.,  NH). 

2.4.  Data  Analysis 

Statistical  analysis  was  carried  out  on  two  types  of  data  inputs:  (a)  the  fluorescence  and  diffuse  reflectance 
spectra  recorded  at  each  of  the  three  illumination-collection  separations,  and  (b)  the  fluorescence  (at  each 
excitation  wavelength)  or  diffuse  reflectance  intensity  integrated  over  the  entire  spectrum  as  a  function  of 
illumination-collection  separation.  A  multivariate  statistical  algorithm  similar  to  the  one  developed  previously 
by  our  group  10  was  used  for  input  (a).  All  data  analyses  were  carried  out  using  a  self-programmed  MATLAB 
code. 

2.4.1.  Spectral  analysis  at  each  of  the  three  illumination-collection  separations 

Briefly,  the  multivariate  statistical  algorithm  consists  of  four  steps  10.  First,  the  spectral  data  are  pre- 
processed  to  minimize  inter-patient  variations  and  variations  due  to  probe-tissue  contact,  by  normalizing  each 
spectrum  to  a  specific  intensity-wavelength  point  in  the  spectra  or  the  integrated  spectral  intensity.  Particularly, 
the  fluorescence  spectra  were  each  normalized  to  their  peak  intensity  and  the  diffuse  reflectance  spectra  were 
normalized  to  their  integrated  intensity.  Second,  the  spectral  data  are  dimensionally  reduced  using  a  multivariate 
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^analysis  technique,  such  as  Principal  Component  Analysis  18.  The  third  step  is  to  perform  feature  extraction,  i.e. 
identify  the  principal  components  (PCs)  that  show  the  statistically  most  significant  differences  between 
malignant  and  non-malignant  breast  tissues  using  a  Wilcoxon  rank-sum  test  19.  The  fourth  step  is  to  use  a 
Support  Vector  Machine  (SVM)  algorithm  20  to  classify  each  sample  as  malignant  or  non-malignant  based  on 
their  PC  scores. 

Data  obtained  in  the  clinical  study  can  be  divided  into  “independent”  and  “dependent”  variables.  The 
independent  variables  are  in  the  spectral  data  matrix  X,  where  each  row  corresponds  to  the  preprocessed 
fluorescence  or  reflectance  spectrum  of  a  tissue  sample  and  each  column  corresponding  to  preprocessed  spectral 
intensity  at  a  specific  wavelength  and  the  “dependent”  variable  Y  is  a  binary  variable  that  represents  the 
histological  diagnosis  of  each  sample,  with  “1”  for  malignant  tissues  and  “0”  for  non-malignant  (benign  and 
normal)  tissues.  Two  multivariate  statistical  analysis  techniques,  Principal  Component  Analysis  (PCA)  18  and 
Partial  Least  Square  (PLS)  regression  21  were  employed  for  spectral  data  reduction.  Both  methods  project  the 
set  of  spectra  onto  a  subspace  expanded  by  the  principal  components  in  order  to  represent  the  spectral  data  with 
a  few  variables,  which  account  for  most  of  the  variance  in  the  original  spectral  data  set. 

In  PCA,  the  set  of  independent  variables  (X)  were  used  for  extraction  of  the  PCs,  that  is,  spectra  measured 
from  all  tissue  samples  regardless  of  histological  diagnoses  were  pooled  together  to  extract  the  PCs.  The  PCs 
were  calculated  such  that  the  first  principal  component  (PCI)  accounts  for  the  largest  amount  of  the  total 
variance  of  the  input  data.  The  second  PC  (PC2)  accounts  for  the  second  largest  amount  of  the  variance  while 
being  orthogonal  to  PCI,  and  so  on. 

PLS  is  a  regression  method  2I,  where  both  the  independent  (X)  and  dependent  (Y)  variables  are  utilized  in 
calculating  the  PCs.  The  PLS  regression  procedure  searches  for  a  set  of  components  that  performs  a 
simultaneous  decomposition  of  X  and  Y  with  the  constraint  that  these  components  explain  as  much  of  the 
covariance  between  X  and  Y.  This  amounts  to  determining  two  sets  of  weights  w  and  c  in  order  to  create  a 
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.linear  aombination  of  the  columns  of  X  and  Y  such  that  their  covariance  is  maximal.  Specifically,  a  pair  of 
vectors  t  =  Xw  and  u  =  Yc  are  obtained  with  the  constraint  that  tTu  is  maximal  (where  the  superscript T  denotes 
the  matrix  transpose).  When  t  has  converged,  it  is  used  to  compute  the  principal  component,  p  for  X,  where  p  = 
XTt  (again  the  superscript  T  denotes  the  matrix  transpose).  When  the  first  PC  (PCI)  is  found,  the  spectral  data 
explained  by  this  PC  is  subtracted  from  both  X  and  Y,  and  the  procedure  is  repeated.  That  is,  additional  PCs  are 
calculated  from  residuals,  which  account  for  the  portion  of  input  data  that  were  not  accounted  for  by  the 
previously  extracted  PCs. 

The  use  of  dependent  variables  (histological  diagnoses)  in  PLS  gives  rise  to  the  difference  between  PCA 
and  PLS.  PCA  provides  PCs  that  account  for  most  of  the  spectral  variance;  however  some  of  these  components 
may  not  necessarily  be  relevant  to  the  diagnoses  of  the  tissue  samples.  PLS  looks  for  PCs  that  not  only  account 
for  a  large  amount  of  the  variance,  but  also  those  which  are  most  relevant  to  the  known  histology  of  the 
samples. 

In  this  study,  PCA  and  PLS  were  performed  on  all  fluorescence  spectra,  one  excitation  wavelength  at  a  time 
and  separately  on  the  diffuse  reflectance  spectra  to  generate  a  set  of  PCs  for  each  spectral  data  set.  Rather  than 
retaining  all  of  the  PCs,  only  a  subset  of  PCs  was  retained  for  further  data  processing.  In  the  case  of  PCA,  the 
PCs  that  account  for  95%  of  the  total  spectral  variance  were  retained  for  further  analysis.  In  the  case  of  PLS,  the 
first  five  PCs  of  each  spectrum  were  retained,  which  account  for  95  -  99%  variance  in  the  spectral  data.  In  both 
cases,  the  original  spectrum  could  be  faithfully  approximated  using  the  linear  combination  of  the  selected  subset 
of  PCs. 

The  scores  of  the  selected  PCs  extracted  from  either  fluorescence  or  diffuse  reflectance  spectra,  were  then 
pooled  together  for  a  Wilcoxon  rank-sum  test.  The  Wilcoxon  rank-sum  test  was  used  to  determine  which  PCs 
showed  the  statistically  most  significant  differences  between  malignant  and  non-malignant  breast  tissues.  The 
scores  of  three  statistically  most  significant  PCs  obtained  using  either  PCA  or  PLS  were  retained  as  inputs  for 
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plassiffcation.  Support  Vector  Machines  (SVM),  including  linear  SVM  and  non-linear  SVM  (e.g.  polynomial 
SVM  and  Radial  Basis  Function  (RBF)  SVM) 20  were  employed  as  classification  schemes  to  specify  a  particular 
sample  as  malignant  or  non-malignant.  Finally,  the  unbiased  performance  of  the  PCA-SVM  and  PLS-SVM 
algorithms  was  evaluated  using  a  “leave  one  out”  scheme  . 

2.4.2.  Integrated  spectral  intensity  analysis 

Integrated  spectral  intensities  were  first  obtained  by  integrating  the  intensities  over  the  entire  spectrum  for 
each  fluorescence  spectrum  and  the  diffuse  reflectance  spectrum  at  a  specific  collection  ring.  Then  intensity 
ratios  were  calculated  by  normalizing  the  integrated  intensity  recorded  from  each  collection  ring  to  that  from  the 
inner  collection  ring  (shortest  illumination-collection  separation).  Thereby  two  intensity  ratios  were  obtained  for 
each  spectrum,  i.e.,  intensity  ratio  of  the  middle  and  inner  ring  (middle-to-inner  intensity  ratio),  and  intensity 
ratio  of  the  outer  and  inner  ring  (outer-to-inner  intensity  ratio).  The  middle-to-inner  and  the  outer-to-inner 
intensity  ratios  of  each  spectrum  (two  variables)  were  then  input  to  the  SVM  classifier  for  discriminating 
between  malignant  and  non-malignant  tissues.  An  unbiased  evaluation  on  the  algorithm’s  performance  was 
obtained  using  the  “leave-one-out”  scheme  . 

3.  Results 

Figure  2  shows  the  average  Excitation  Emission  Matrices  (EEMs)  of  malignant  (n  =  13),  normal/benign 

fibrous  (n  =  14)  and  normal  adipose  (n  =  20)  tissues.  Figures  (a)  -  (c)  correspond  to  the  average  EEMs  of 

malignant  tissues  measured  with  the  inner  (a),  middle  (b),  and  outer  (c)  collection  rings,  respectively,  and  figure 

(d)  and  (e)  correspond  respectively  to  average  EEMs  of  normal/benign  fibrous  and  normal  adipose  tissues 

measured  with  the  inner  collection  ring.  All  figures  are  plotted  on  a  log  scale.  Average  EEMs  of  malignant 

tissues  (Figs.  2  (a)-(c))  show  that  four  peaks  are  visible  and  appear  at  similar  locations  in  all  three  EEMs,  i.e.  at 

excitation-emission  wavelength  pairs  of  300-340,  340-390,  360/380-460  and  440-520  nm.  As  expected,  the 

fluorescence  intensity  at  each  peak  location  decreases  from  the  inner  ring  to  the  outer  ring.  Such  a  decrease  in 

fluorescence  intensity  is  also  observed  in  the  EEMs  of  normal/benign  fibrous  and  normal  adipose  tissues 
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rpeasur&d  with  three  collection  rings  (not  shown  here).  The  average  EEM  of  normal/benign  fibrous  tissues 
measured  with  the  inner  ring  (Fig.  2  (d))  has  the  same  four  peaks  as  the  malignant  tissue  EEM  measured  with 
the  inner  ring  (Fig.  2  (a)).  However,  the  average  adipose  tissue  EEM  measured  with  the  inner  ring  (Fig.  2  (e)) 
displays  distinct  differences,  with  a  weak  presence  of  the  340-390  nm  peak  and  a  red  shift  in  the  peak  at  360- 
460  nm  to  approximately  360-500  nm. 


(Figure  2) 

Figure  3  displays  (a)  the  average  diffuse  reflectance  spectra  of  malignant  breast  tissues  (n  =  13)  measured 
with  the  inner,  middle  and  outer  rings,  and  (b)  the  average  diffuse  reflectance  spectra  of  malignant  (n  =  13), 
normal/benign  fibrous  (n  =  14)  and  normal  adipose  tissues  (n  =  20)  measured  with  the  inner  ring.  In  each 
spectrum,  five  valleys  are  notable  at  350,  420,  470,  540  and  570  nm.  The  valleys  at  420,  540  and  570  nm 
correspond  to  hemoglobin  absorption,  the  one  at  350  nm  corresponds  to  NADH  and  the  one  at  470  nm  can  be 
assigned  to  P-carotene.  In  Fig.  3(a),  the  diffuse  reflectance  intensity  measured  from  malignant  breast  tissues 
decreases  from  the  inner  collection  ring  to  the  outer  collection  ring,  particularly  in  the  wavelength  range  above 
450  nm.  Such  a  decrease  is  also  observed  in  diffuse  reflectance  spectra  of  normal/benign  fibrous  and  normal 
adipose  tissues  measured  with  the  three  collection  rings  (not  shown  here).  A  comparison  of  spectral  intensities 
for  the  three  tissue  types,  measured  with  the  inner  ring  (Fig.  3(b))  shows  that  fibrous  tissues  have  the  highest 
average  diffuse  reflectance  intensity,  while  adipose  tissues  have  the  lowest  average  diffuse  reflectance  intensity. 

(Figure  3) 

PCA  and  PLS  were  carried  out  on  the  fluorescence  and  diffuse  reflectance  spectra  measured  with  each  of  the 
three  collection  rings,  to  obtain  a  set  of  PCs  that  faithfully  represent  the  spectral  data.  Table  3  shows  the  three 
statistically  most  significant  PCs  for  each  collection  ring,  identified  by  a  Wilcoxon  rank-sum  test  from  the  two 
sets  of  PCs  obtained  from  the  PCA  and  PLS  analyses.  All  the  PCs  here  display  statistically  significant  differences 
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between  malignant  and  non-malignant  tissues  below  a  significance  level  of  p  <  0.005.  The  three  most  significant 
PCs  obtained  from  PLS  analysis  were  extracted  from  fluorescence  spectra  at  320  and  420  nm  excitation  and  from 
the  diffuse  reflectance  spectra.  The  three  most  significant  PCs  obtained  from  PCA  analysis  were  extracted  from 
fluorescence  spectra  at  320,  420,  300  and  340  nm  excitation  and  from  the  diffuse  reflectance  spectra. 

(Table  3) 

Table  4  shows  the  overall  classification  rate,  sensitivity  and  specificity  achieved  with  linear  SVM  for 
discriminating  between  malignant  and  non-malignant  breast  tissues  using  the  three  statistically  most  significant 
PCs  (obtained  from  PCA  and  PLS  analysis)  for  each  collection  ring.  It  should  be  noted  that  in  several  cases, 
particularly  in  the  case  where  PLS  analysis  was  used,  the  three  most  significant  PCs  included  those  obtained  from 
both  diffuse  reflectance  and  fluorescence  spectra.  The  classification  accuracy  achieved  with  PCA-SVM  or  PLS- 
SVM  did  not  differ  significantly  between  the  three  collection  rings.  However,  for  spectral  data  acquired  by  a 
given  collection  ring,  the  PLS-SVM  algorithm  achieved  consistently  higher  classification  accuracy  than  the  PCA- 
SVM  algorithm.  It  was  also  found  that  the  overall  classification  rate,  sensitivity  and  specificity  obtained  using 
non-linear  SVM  classification  (polynomial  SVM  and  Radial  Basis  Function  (RBF)  SVM)  did  not  differ 
significantly  from  that  obtained  using  linear  SVM  classification  (not  shown  here). 

(Table  4) 

Table  5  shows  the  overall  classification  rate,  sensitivity  and  specificity  achieved  with  linear  SVM  for 

discriminating  between  malignant  and  non-malignant  breast  tissues  using  the  three  statistically  most  significant 

«’ 

PCs  extracted  from  PLS  analysis  of  the  diffuse  reflectance  spectra  measured  with  each  collection  ring.  The  three 
PCs  show  statistically  significant  differences  between  malignant  and  non-malignant  tissues  below  a  significance 
level  of  p  <  0.05.  Only  the  results  from  the  PLS-SVM  algorithm  are  shown,  since  this  approach  yields  superior 
results  compared  to  the  PCA-SVM  algorithm.  The  overall  classification  rate  did  not  differ  significantly  between 
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.the  different  collection  rings.  In  addition,  the  specificity  was  higher  than  the  sensitivity  for  all  three  collection 
rings.  A  comparison  of  the  classification  accuracy  achieved  using  (a)  combined  fluorescence  and  reflectance 
spectra  (Table  4),  versus  (b)  reflectance  spectra  only  (Table  5)  indicates  that  the  overall  classification  rates  do  not 
differ  significantly  for  these  two  cases.  However,  the  sensitivity  is  consistently  higher  for  the  former,  while  the 
specificity  is  consistently  higher  for  the  latter. 


(Table  5) 

Figure  4  shows  the  average  integrated  spectral  intensity  ratios  and  standard  deviations  for  (a)  fluorescence 
spectra  at  300  nm,  excitation,  and  (b)  diffuse  reflectance  spectra  of  malignant  (n  =  13),  normal/benign  fibrous  (n 
=  14)  and  normal  adipose  tissues  (n  =  20).  For  fluorescence  at  300  nm  excitation,  the  average  values  of  middle-to- 
inner  ratio  and  outer-to-inner  ratio  were  0.29  +  0.12  and  0.21  ±  0.06  respectively  for  malignant  tissues,  0.35  ± 
0.08  and  0.25  ±  0.04  for  normal/benign  fibrous  tissues,  and  0.49  ±  0.09  and  0.35  ±  0.06  for  normal  adipose 
tissues.  For  the  diffuse  reflectance,  the  average  values  of  middle-to-inner  ratio  and  outer-to-inner  ratio  were  0.52 
±  0.08  and  0.48  ±  0.07  respectively  for  malignant  tissues,  0.48  ±  0.07  and  0.43  ±  0.05  for  normal/benign  fibrous 
tissues,  and  0.62  ±  0.08  and  0.59  ±0.12  for  normal  adipose  tissues.  As  expected,  the  spectral  intensity  ratios  of 
both  fluorescence  and  diffuse  reflectance  decrease  with  increasing  illumination-collection  separation. 
Furthermore  the  decrease  (attenuation)  in  spectral  intensity  is  greater  from  inner  ring  to  middle  ring  than  from 
middle  ring  to  outer  ring.  Also,  the  attenuation  in  spectral  intensity  is  different  for  different  tissue  types.  In 
particular,  normal  adipose  tissues  undergo  less  attenuation  relative  to  normal/benign  fibrous  and  malignant 
tissues.  Spectral  intensity  ratios  at  the  other  seven  excitation  wavelengths  are  not  shown  here,  but  display  a 
similar  pattern  to  that  observed  at  an  excitation  wavelength  of  300  nm.  Wilcoxon  rank-sum  tests  showed  that  the 
middle-to-inner  and  outer-to-inner  intensity  ratios  of  fluorescence  at  300  nm,  excitation,  displayed  statistically 
significant  differences  between  malignant  and  non-malignant  tissues  below  a  significance  level  of  p  <  0.005.  The 
same  test  indicated  that  the  middle-to-inner  and  outer-to-inner  intensity  ratios  of  diffuse  reflectance  did  not 

display  statistically  significant  differences  between  malignant  and  non-malignant  tissues  (p  <  0.5). 
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(Figure  4) 


Table  6  shows  the  overall  classification  rate,  sensitivity  and  specificity  achieved  with  linear  SVM  for 
discriminating  between  malignant  and  non-malignant  tissues  using  the  middle-to-inner  and  outer-to-inner 
intensity  ratios  for  fluorescence  collected  at  different  excitation  wavelengths  and  for  the  diffuse  reflectance.  The 
results  show  that  classification  using  fluorescence  spectral  intensity  ratios  at  300  nm  excitation  provides  the 
highest  classification  accuracy,  which  is  comparable  to  that  obtained  using  the  entire  spectral  input  for  a  given 
collection  ring  (Table  4).  Using  fluorescence  at  320  to  360  nm  excitation  provides  slightly  lower  classification 
accuracy,  which  is  still  comparable  to  that  obtained  using  the  entire  spectral  input  for  a  given  collection  ring. 
When  spectral  intensity  ratios  of  either  fluorescence  at  excitation  wavelengths  of  380  to  440  nm,  or  diffuse 
reflectance  were  used,  the  classification  yields  high  specificity,  but  very  low  sensitivity. 

(Table  6) 

Figure  5  shows  the  scatter  plots  of  the  middle-to-inner  and  outer-to-inner  ratios  for  (a)  fluorescence  spectra 
at  300  nm  excitation  and  (b)  diffuse  reflectance  spectra  of  malignant  and  non-malignant  breast  tissues  and  the 
corresponding  hyperplanes  obtained  from  linear  SVM.  In  the  case  of  fluorescence  at  300  nm  excitation  (Fig. 
5(a)),  the  clusters  of  the  spectral  intensity  ratios  for  the  two  tissue  types  are  separable.  However  in  the  case  of 
diffuse  reflectance  (Fig.  5(b),  the  clusters  of  spectral  intensity  ratios  for  the  two  tissue  types  are  non-separable. 
Similar  observations  were  made  from  scatter  plots  of  spectral  intensity  ratios  for  fluorescence  spectra  at 
excitation  wavelengths  of  380  to  440  nm  (not  shown  here). 

4.  Discussion 

Optical  spectroscopy  in  the  UV-VIS  spectral  range  shows  distinct  differences  between  malignant  and  non- 

malignant  breast  tissues,  and  thus  this  technique  has  the  potential  to  be  used  as  a  diagnostic  tool  for  breast 
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.cancer,.  Statistical  analysis  of  the  spectra  indicated  that  the  fluorescence  spectra  at  excitation  wavelengths  of 
300,  320,  340  and  420  nm  and  the  diffuse  reflectance  spectrum  showed  the  statistically  most  significant 
differences  (p  <  0.005)  between  malignant  and  non-malignant  breast  tissues.  A  previous  study  by  our  group  10 
showed  that  fluorescence  spectra  at  excitation  wavelengths  of  300,  320,  400,  420,  440  and  460  nm  showed  the 
statistically  most  significant  differences  between  malignant  and  non-malignant  breast  tissues  (p  <  0.005).  The 
excitation  wavelengths  identified  in  the  current  study  are  a  subset  of  the  wavelengths  identified  previously, 
except  for  the  excitation  wavelength  of  340  nm.  Another  difference  is  that  in  the  previous  study,  the  diffuse 
reflectance  spectrum  was  not  found  to  show  statistically  significant  differences  between  malignant  and  non- 
malignant  breast  tissues. 

Two  statistical  methods,  PCA  and  PLS,  were  used  in  this  study  for  spectral  data  reduction.  One  of  the 
differences  between  PCA  and  PLS  is  that  in  PCA,  the  components  are  chosen  so  that  maximal  data  variance  is 
explained,  while  PLS  looks  for  components  that  not  only  describe  as  much  of  the  data  variance  as  possible,  but 
those  which  are  also  most  relevant  to  the  known  group  association  of  each  sample.  A  linear  SVM  algorithm 
based  on  PCs  obtained  from  PLS  outperformed  that  based  on  PCs  obtained  from  PCA  for  spectra  collected  at  all 
three  illumination-collection  separations  (Table  4).  In  order  to  determine  how  different  the  PCs  obtained  from 
PCA  and  PLS  are  for  a  given  spectrum  measured  at  a  particular  illumination-collection  separation,  those 
extracted  from  fluorescence  spectra  at  320  and  420  nm  excitation  and  from  the  diffuse  reflectance  spectra 
measured  with  the  inner  collection  ring  of  the  fiber-optic  probe  (see  first  row  of  Table  3)  were  compared. 
Although  the  line  shapes  of  the  PLS  and  PCA  components  extracted  from  the  fluorescence  and  diffuse 
reflectance  spectra  were  found  to  be  generally  similar,  the  PLS  components  had  additional  structural  features 
that  were  not  present  in  the  PCA  components.  However,  it  is  difficult  to  directly  assign  these  features  to 
particular  fluorophores  and  chromophores  in  the  tissue.  This  is  due  to  the  fact  that  these  linear  models  do  not 
accurately  describe  the  non-linear  relationship  between  fluorescence,  absorption  and  scattering  in  turbid  media 
such  as  tissue.  Moreover,  a  component  can  be  multiplied  by  -1  without  affecting  the  decomposition  model. 
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Thus,  peaks  in  the  PCA/PLS  components  could  correspond  to  peaks  or  valleys  in  a  given  spectrum,  making  it 
difficult  to  assign  the  spectral  features  to  particular  fluorophores  or  chromophores  in  the  tissue. 

Typically,  the  optical  spectral  features  of  different  tissue  types  are  not  completely  separable,  and  a  non¬ 
linear  SVM  is  expected  to  be  superior  to  a  linear  SVM  classification  in  dealing  with  such  non-separable  cases. 
However,  using  the  same  data  reduction  technique,  PCA  or  PLS,  the  classification  rates,  sensitivities  and 
specificities  provided  by  linear  and  non-linear  SVMs  did  not  differ  significantly.  On  the  other  hand,  linear  SVM 
classification  associated  with  PLS  data  reduction  (PLS-SVM)  provided  higher  classification  accuracy  than 
linear  SVM  associated  with  PCA  data  reduction  (PCA-SVM)  (Table  4).  This  suggests  that  the  improvement  in 
the  classification  accuracy  is  more  heavily  dependent  on  the  data  reduction  technique.  An  alternative  approach 
to  improving  the  classification  accuracy  could  be  to  employ  physical  models  to  extract  tissue  optical  properties 
that  are  diagnostically  useful  23 .  This  will  be  explored  in  future  studies. 

In  this  study,  fluorescence  and  diffuse  reflectance  spectra  were  measured  from  breast  tissues  with  a  multi¬ 
separation  fiber-optic  probe.  This  probe  has  three  illumination-  collection  separations,  which  enable  optical 
spectra  to  be  measured  from  different  depths  within  the  tissue.  The  classification  accuracy  achieved  with  PCA- 
SVM  or  PLS-SVM  did  not  differ  significantly  between  the  three  illumination-collection  separations.  This  is 
likely  due  to  the  fact  that  (1)  the  probing  depth  is  expected  to  increase  only  modestly  with  increasing 
illumination-collection  separation  (based  on  the  Monte  Carlo  simulation  results  in  Table  2)  and  (2)  there  is  no 
apparent  structural  change  in  breast  tissue  with  increasing  depth  that  could  be  captured  with  different 
illumination-collection  separations. 

The  fluorescence/reflectance  signal  attenuates  with  increasing  probing  depth  in  tissue  (due  to  the  increased 
path  length  of  the  photons)  and  the  multi-separation  fiber-optic  probe  can  be  used  to  measure  the  attenuation 
characteristics  of  these  signals  in  the  breast.  The  results  from  the  data  analysis  show  that  using  the  integrated 
fluorescence  intensities  recorded  at  a  single  excitation  wavelength  at  all  three  illumination-collection 
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^separations  (a  measure  of  attenuation)  can  discriminate  malignant  from  non-malignant  breast  tissues  with 
similar  classification  accuracy  to  that  using  spectral  data  measured  at  several  excitation  wavelengths  at  a  single 
illumination-collection  separation.  The  use  of  integrated  fluorescence  intensities  at  a  single  excitation 
wavelength  is  expected  to  (1)  significantly  reduce  the  measurement  time  and  (2)  require  much  less  data 
processing.  This  finding  has  significant  implications  in  clinical  applications,  where  both  speed  and  low  cost  are 
desirable.  Although  each  spectrum  takes  only  a  few  seconds  to  measure,  hardware  changes  between 
measurements  increases  the  time  for  a  single  scan  of  8  fluorescence  spectra  and  one  diffuse  reflectance 
spectrum  to  a  minute.  If  the  system  can  be  streamlined  to  include  a  subset  of  the  wavelengths  evaluated  in  this 
study  (such  as  the  integrated  intensity  at  a  single  excitation  wavelength),  the  measurement  and  data  processing 
time  can  be  significantly  reduced  to  enable  near  real-time  implementation  of  this  technology  during  a  clinical 
procedure. 

Note  that  the  number  of  non-malignant  samples  (n  =  34)  is  much  greater  than  that  of  malignant  samples  (n  = 
13)  in  this  study.  The  principal  idea  of  an  SVM  is  to  determine  an  optimal  hyperplane  that  maximizes  the 
margin  between  two  classes,  and  equalize  the  distances  of  the  misclassified  samples  in  each  class  from  the 
separation  boundary.  For  a  poorly  separable  case  (as  in  Fig.  5(b))  and  when  the  numbers  of  elements  in  the  two 
classes  are  unbalanced,  a  large  amount  of  samples  from  each  class  become  support  vectors  (i.e.  marginal 
samples),  and  it  is  likely  that  most  samples  from  the  class  which  has  much  fewer  elements  (in  this  case,  the 
malignant  tissues)  will  be  misclassified.  This  would  result  in  a  high  specificity,  but  very  low  sensitivity. 
However  for  the  case  where  sample  clusters  are  separable  (Fig.  5  (a)),  it  is  less  likely  that  the  sensitivity  and 
specificity  are  sensitive  to  the  number  of  samples  in  each  class.  This  also  likely  explains  the  high  specificity  and 
low  sensitivity  achieved  when  diffuse  reflectance  spectra  alone  (Table  5)  rather  than  the  combination  of 
fluorescence  and  diffuse  reflectance  spectra  (Table  4)  were  used  for  discriminating  between  malignant  and  non- 
malignant  breast  tissues.  The  three  PCs  extracted  from  diffuse  reflectance  spectra  showed  statistically  smaller 
differences  between  malignant  and  non-malignant  tissues  (p  <  0.05)  compared  to  the  three  PCs  extracted  from 
fluorescence  and  diffuse  reflectance  spectra  (p  <  0.005).  To  verify  this  explanation,  data  analyses  (using  spectral 
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intensity  ratios  only)  were  carried  out  on  a  subset  of  spectral  data,  which  included  balanced  malignant  and  non- 
malignant  samples.  In  this  data  set,  only  half  of  the  normal/benign  fibrous  samples  (7)  and  half  of  normal 
adipose  samples  (10)  were  retained,  and  the  number  of  malignant  samples  vs.  non-malignant  samples  was 
13:17.  The  classification  results  indicated  that  for  the  separable  case  (i.e.  fluorescence  at  an  excitation 
wavelength  of  300  nm),  the  sensitivity  and  specificity  achieved  with  the  balanced  data  set  (61.54%  and  82.35% 
respectively)  are  comparable  to  that  achieved  with  the  unbalanced  data  set  (69.23%  and  88.24%  respectively). 
However  for  the  non-separable  case,  there  was  a  significant  disparity  in  the  sensitivities  and  specificities 
achieved  with  the  balanced  and  unbalanced  data  sets  (53.85%  and  68.75%  for  the  balanced  data  set  vs.  0%  and 
100%  for  unbalanced  data  set).  In  summary,  when  sample  clusters  are  separable,  the  sensitivity  and  specificity 
do  not  appear  to  be  highly  sensitive  to  the  number  of  samples  in  each  class.  However  when  samples  are  poorly 
separable,  it  is  likely  the  specificity  overwhelms  the  sensitivity  when  there  are  too  many  non-malignant  vs. 
malignant  samples. 

The  collection  of  tissue  samples  in  this  study  provided  a  relatively  small  data  set,  with  a  total  of  47  samples, 
only  13  of  which  were  malignant.  Although  a  “leave-one-out”  cross  validation  scheme  provides  an  unbiased 
evaluation  of  the  algorithm,  a  larger  sample  size  will  need  to  be  accrued  in  future  studies  to  fully  validate 
classification  results  reported  in  this  paper.  However,  it  should  be  noted  that  similar  classification  accuracies 
are  obtained  in  this  study  when  different  preprocessing  techniques  (spectra  vs.  integrated  spectral  intensity), 
feature  extraction  methods  and  classifiers  (linear  and  non-linear  SVM)  are  used.  These  results  suggest  that  the 
difference  in  the  spectral  features  between  malignant  and  non-malignant  tissues  can  be  identified  consistently 
using  a  variety  of  different  algorithms. 

The  study  presented  here  and  our  previous  study  10  provide  the  foundation  for  the  primary  focus  of  the 
future  work,  which  is  to  explore  the  utility  of  optical  spectroscopy  for  breast  cancer  detection  during  core  needle 
biopsy.  The  current  probe  geometry  will  be  modified  into  a  side-firing  fiber-optic  probe  and  implemented  for 
use  in  a  vacuum-assisted  core  biopsy  needle.  If  this  optical  technique  proves  to  be  diagnostically  useful,  it  can 
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potentially  improve  the  diagnostic  efficacy  of  breast  needle  biopsy,  and  lead  to  fewer  biopsies  and  follow  up 
procedures  in  patients  suspected  to  have  breast  cancer.  Additionally,  the  fiber-optic  probe  can  be  made  thin 
enough  to  fit  through  an  even  smaller  needle  than  the  standard  1 /4-inch  size,  making  an  emotionally  draining 
procedure  less  physically  traumatic. 
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Table  Captions 

Table  1.  The  histological  breakdown  of  the  47  samples  examined  in  the  breast  tissue  optical  spectroscopy 
study. 

Table  2.  The  80%  probing  depth  and  the  percent  fluorescence  detected  with  the  three  illumination-collection 
pairs  of  the  multi-separation  probe  from  a  homogeneous  fluorescent  medium,  with  (a)  variable  absorption 
coefficients  (pa)  and  a  fixed  scattering  coefficient  (ps),  and  (b)  variable  scattering  coefficients  and  a  fixed 
absorption  coefficient.  The  anisotropy  factor  (g)  of  the  medium  was  set  to  0.926.  Note:  QY  is  quantum  yield, 
Xexc  is  the  excitation  wavelength  and  Xemm  is  the  emission  wavelength. 

Table  3.  The  three  statistically  most  significant  PCs  for  each  collection  ring,  identified  by  a  Wilcoxon  rank-sum 
test  from  the  two  sets  of  PCs  obtained  from  the  PCA  and  PLS  analyses.  All  the  PCs  here  display  statistically 
significant  differences  between  malignant  and  non-malignant  tissues  below  a  significance  level  of  p  <  0.005. 

Table  4.  The  overall  classification  rate,  sensitivity  and  specificity  achieved  with  linear  SVM  for  discriminating 
between  malignant  and  non-malignant  breast  tissues  using  the  three  statistically  most  significant  PCs  (obtained 
from  PCA  and  PLS  analysis)  for  each  collection  ring. 

Table  5.  The  overall  classification  rate,  sensitivity  and  specificity  achieved  with  linear  SVM  for  discriminating 
between  malignant  and  non-malignant  breast  tissues  using  the  three  statistically  most  significant  PCs  extracted 
from  PLS  analysis  of  the  diffuse  reflectance  spectra  measured  with  each  collection  ring. 

Table  6.  The  overall  classification  rate,  sensitivity  and  specificity  achieved  with  linear  SVM  for  discriminating 
between  malignant  and  non-malignant  tissues  using  the  middle-to-inner  and  outer-to-inner  intensity  ratios  for 
fluorescence  collected  at  different  excitation  wavelengths  and  for  the  diffuse  reflectance. 
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Figure  Captions 


Figure  1.  (a)  Schematic  of  the  optical  spectrometer  and  (b)  common  end  of  the  fiber-optic  probe  that  comes  in 
contact  with  the  tissue  (the  black  circles  correspond  to  the  illumination  fibers,  the  gray  circles  correspond  to  the 
collections  fibers  which  form  the  three  collection  rings,  and  the  white  circles  correspond  to  dead  fibers  for 
bundle  packing). 

Figure  2.  The  average  Excitation  Emission  Matrices  (EEMs)  of  malignant  (n  =  13),  normal/benign  fibrous  (n  = 
14)  and  normal  adipose  (n  =  20)  tissues.  Figures  (a)  -  (c)  correspond  to  the  average  EEMs  of  malignant  tissues 
measured  with  the  inner  (a),  middle  (b),  and  outer  (c)  collection  rings,  respectively  and  figure  (d)  and  (e) 
correspond  respectively  to  average  EEMs  of  normal/benign  fibrous  and  normal  adipose  tissues  measured  with 
the  inner  collection  ring.  All  figures  are  plotted  on  a  log  scale. 

Figure  3.  (a)  The  average  diffuse  reflectance  spectra  of  malignant  breast  tissues  (n  =  13)  measured  with  the 
inner,  middle  and  outer  rings,  and  (b)  the  average  diffuse  reflectance  spectra  of  malignant  (n  =  13), 
normal/benign  fibrous  (n  =  14)  and  normal  adipose  tissues  (n  =  20)  measured  with  the  inner  ring. 

Figure  4.  Average  integrated  spectral  intensity  ratios  and  standard  deviations  for  (a)  fluorescence  spectra  at  300 
nm,  excitation,  and  (b)  diffuse  reflectance  spectra  of  malignant  (n  =  13),  normal/benign  fibrous  (n  =  14)  and 
normal  adipose  tissues  (n  =  20). 

Figure  5.  Scatter  plots  of  middle-to-inner  and  outer-to-inner  intensity  ratios  for  (a)  fluorescence  spectra  at  300 
nm,  excitation,  and  (b)  diffuse  reflectance  spectra  of  malignant  and  non-malignant  breast  tissues  and  the 
corresponding  hyperplanes  obtained  from  linear  SVM. 
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Tables* 


Table  1. 


Histological  Breakdown 

Subcategory 

Number  of  Samples 

Invasive  ductal  carcinoma  (IDC) 

7 

Malignant  tissues 

Invasive  lobular  carcinoma  (ILC) 

4 

13 

Ductal  carcinoma  in  situ  (DCIS) 

2 

Normal 

Fibrous 

12 

32 

Adipose 

20 

Benign 

2 

2 

Total 

47 
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Table  2# 


& 


Optical  properties 

80%  probing  depth  (pm) 

Percent  fluorescence 

^-exc 

ps  =  110.4  cm-1 

(x  10'4) 

Pa 

(cm'1) 

QY 

H 

Inner 

Middle 

Outer 

Inner 

Middle 

Outer 

1.312 

0.762 

0.994 

HI 

■a 

IBB 

14.7 

10.1 

7.62 

2.025 

0.494 

1.596 

in 

■EES 

ns 

17.7 

11.6 

8.46 

5.286 

0.189 

4.352 

900 

1000 

1100 

21.6 

11.7 

7.16 

10.316 

0.097 

8.580 

700 

800 

900 

18.5 

7.92 

4.00 

14.427 

0.069 

12.012 

650 

700 

800 

15.4 

5.28 

2.30 

19.191 

0.052 

16.043 

550 

600 

700 

11.9 

3.54 

1.26 

21.236 

0.047 

17.754 

550 

550 

550 

10.6 

2.72 

1.03 

31.812 

0.031 

26.517 

400 

500 

500 

6.38 

1.23 

0.292 

lb) 


Optical  properties 

80%  probing  depth  (pm) 

Percent  fluorescence 

^■cxc 

pa=  10.82  cm"1 

^emm 

pa=  9.0  cm'1 

(x  10'4) 

Ps 

_ (cm'  ) _ 

QY 

HI 

Inner 

Middle 

Outer 

Inner 

Middle 

Outer 

1 

47.2 

mm 

1000 

12.3 

4.60 

2.34 

1 

70.8 

■n 

15.7 

5.90 

2.80 

100.0 

1 

94.4 

700 

750 

900 

17.0 

6.72 

3.50 

125.0 

1 

118.0 

650 

700 

800 

19.7 

8.06 

4.00 

150.0 

1 

141.6 

600 

650 

800 

21.4 

8.70 

4.14 

175.0 

1 

165.3 

600 

600 

700 

22.2 

9.66 

4.62 

1 

188.9 

500 

550 

600 

24.4 

9.94 

5.02 

1 

212.5 

500 

550 

550 

25.2 

10.3 

4.90 
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Tabled. 


PCA 

— >  increasing  signi 

Icance 

PLS 

— >  increasing  signii 

"icance 

Inner 

420  nm 
PCI 

320  nm 
PCI 

reflectance 

PC4 

Inner 

420  nm 
PCI 

reflectance 

PCI 

320  nm 
PCI 

Middle 

320  nm 
PCI 

420  nm 
PCI 

340  nm 
PCI 

Middle 

420  nm 
PCI 

reflectance 

PCI 

320  nm 
PCI 

Outer 

320  nm 
PCI 

300  nm 
PC3 

420  nm 
PCI 

Outer 

420  nm 
PCI 

reflectance 

PCI 

320  nm 
PCI 

Table  4. 


PCA  +  Linear  SVM 

PLS  +  Linear  SVM 

Inner 

Middle 

Outer 

Inner 

Middle 

Outer 

Classification  Rate  (%) 

78.72 

65.96 

57.45 

82.98 

82.98 

80.85 

Sensitivity  (%) 

38.46 

30.77 

15.38 

61.54 

69.23 

69.23 

Specificity  (%) 

94.12 

79.41 

73.53 

91.18 

88.24 

85.29 

Table  5. 


PLS  +  Linear  SVM 

Inner 

Middle 

Outer 

Classification  Rate  (%) 

80.85 

78.72 

85.11 

Sensitivity  (%) 

53.85 

46.15 

Specificity  (%) 

91.18 

100 

Table  6. 


300  nm 

320  nm 

340  nm 

360  nm 

380  nm 

400  nm 

420  nm 

440  nm 

Refl 

Classification 
Rate  (%) 

82.98 

78.72 

78.72 

78.72 

68.09 

68.09 

65.96 

61.70 

72.34 

Sensitivity  (%) 

69.23 

61.54 

53.85 

61.54 

23.08 

23.08 

0 

0 

0 

Specificity  (%) 

88.24 

85.29 

88.24 

85.29 

85.29 

85.29 

91.18 

85.29 

100 
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Figures 


Figure  1. 
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Reflectance  Intensity  Ratio  Fluorescence  intensity  Ratio 
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